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An analytical model has been developed to evaluate "displacer gap losses" in the clearance between the 
displacer and the cylinder in both (1 and y configurations of Stirling engines. Displacer gap losses are the 
sum of the "shuttle heat transfer" and the "enthalpy pumping". 

The present model takes into account the pressure gradient in the gap, the gas compressibility and real 
gas effect. Gas velocity and temperature distributions in the displacer, in the gap and in the cylinder were 
determined by solving momentum and energy balances in a concentric tubes geometry. Our model is 
then introduced in a whole engine model and the effect of the clearance thickness and engine's speed on 
total displacer gap losses are investigated. Novel tendencies of the solution are observed and new ways 
for optimization are demonstrated. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

The Stirling engine has been invented since the 19th century but 
was unable to penetrate the energy and transport sector despite its 
flexibility and relatively high thermal efficiency. The use of this 
technology has been concentrated in cryogenics and spatial appli¬ 
cations. However, with the present growth of concerns about en¬ 
ergy saving and environmental issues, Stirling engines becomes 
more interesting than ever. It is the best solution to make micro¬ 
scale co-generation since the ORC (organic Rankine cycle) is rela¬ 
tively less efficient for a power level below 100 kW e i [1], Stirling 
cycle implies that a constant mass of working gas (usually air, he¬ 
lium or hydrogen) alternates between two temperature levels and 
passes through four thermodynamic processes: an isothermal 
compression, an isochoric heating, an isothermal expansion and 
finally an isochoric cooling [2]. In /? and y configurations of Stirling 
engines (Fig. 1) this reciprocating motion is assured by a displacer. 
To avoid frictional losses and excessive wear, clearance seals are 
generally used instead of ring seals. The inconvenience of this so¬ 
lution is that a gas leakage occurs from the compression to the 
expansion chamber or back depending on the displacer's direction 
of motion and the pressure difference between these spaces. This 
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leakage causes an enthalpy pumping through the displacer's 
clearance [3], 

Besides, a temperature gradient exists along the cylinder wall. 
Thus, the displacer is moving between a hot and a cold region and 
causes heat transfer with its motion. It takes some heat from the 
hot region of the cylinder and rejects it in the cold region. This 
amount of lost heat is called "shuttle heat transfer". The sum of 
enthalpy pumping and shuttle heat transfer loss are known in the 
litterature as "Dispalcer gap losses" [3]. 

Urieli [4], used a steady state model (steady displacer velocity 
and axial pressure gradient) to calculate the leakage mass flow 
through the displacer's clearance. Urieli's formula was used for 
several years. Rios [5] developed an approximate solution by line¬ 
arization and application of Fourier series to calculate the shuttle 
heat transfer. Later, Baik and Chang [6] developed an analytic so¬ 
lution for shuttle heat transfer by considering only conduction in 
both cylinder and displacer and by neglecting the effect of the gap 
flow between them. This effect was introduced later by Chang et al. 
[7], But, they considered a simple geometry (parallel plates) and 
assumed that the gas flow in the gap is only due to the displace 
motion without any effect of the pressure gradient in the gap. 
Finally, Kotsubo and Swift [8] applied the thermoacoustic theory to 
study enthalpy pumping and shuttle heat transfer. They considered 
a parallel plates geometry and supposed that the leakage flow is a 
priori known. 

Some numerical studies exist also like the ID study proposed by 
Andersen [9] and the work of Huang and Berggren [10], These 
studies have been conducted for specific engines and used a 
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Nomenclature 

A 

thermal conductivity, W K _1 m _1 



V 

kinematic viscosity, m 2 s _1 

c P 

specific heat capacity, J K _1 kg 1 

(j) 

angular speed, rad°s _1 

L 

length, m 



P 

pressure, Pa 

Subscripts 

Pr 

Prandtl number 

0 

time-average 

R 

radius, m or perfect gas spec const., J kg -1 I< 1 

c 

cylinder or compression space 

r 

radial coordinate, m 

d 

displacer 

T 

temperature, I< 

e 

expansion space 

t 

time, s 

f 

fluid 

u 

velocity, ms 1 

H 

high temperature source 

V 

volume, m 3 

h 

hot HEX 

X 

position, m 

k 

cold HEX 

z 

axial coordinate, m 

L 

low temperature sink 



n 

nth order of Fourier series 

Greek letters 



a 

thermal diffusivity, m 2 s _1 

Superscripts 

p 

isochoric thermal pressure coefficient, K 1 

* 

complex number 

r 

temperature gradient in the gap, Km -1 

~ 

oscillating part 

T 

ratio Cp/C v 




convective heat transfer coefficient calculated from correlations. 
Besides, numeric studies are greedy in terms of computing time 
and then can not be used in whole engine optimization approaches. 

In the present study, efforts are put in modeling analytically 
displacer gap losses by solving momentum and energy equations in 
the annular gap between the cylinder and the displacer. Unlike Baik 
[6] and Chang [7], flow and pressure gradient in the gap are taken 
into account here and that affects considerably solutions and new 
phenomena are observed. Introducing our model in a whole engine 
model, we demonstrate the existence of an optimal clearance 
thickness and an optimal engine speed that minimizes displacer 
gap losses. 



Fig. 1 . (3 and y configurations of the Stirling engine: 1 - expansion chamber; 2 — 
displacer; 3 - heater; 4 - regenerator; 5 - cooler; 6 - compression chamber; 7 — 
power piston; 8 — cylinder. 


2. Problem formulation 

2.1. Problem description 

As shown in Fig. 1, The displacer separates the compression and 
the expansion spaces. The pressure difference between these 
spaces is caused by the friction of the gas flowing through the heat 
exchangers and the regenerator. The problem is modeled in this 
work as a pulsating flow in the annular gap between two concentric 
tubes: the displacer and the cylinder. The displacer is moving 
periodically relative to the cylinder. In Fig. 2, two types of coordi¬ 
nate systems are used: a stationnary system (z,r) for the cylinder 
and the gap fluid and a moving system (z d ,r d ) for the displacer. 
These coordinate systems are related by the displacer’s motion: 

r d = r (la) 

z d =z-x d (t) (lb) 

where x d (t) is the axial coordinate of the origin of the moving co¬ 
ordinate system in the stationary one. It represents also the posi¬ 
tion of the displacer relative to the cylinder. 

2.2. Assumptions 

As said previously, we keep some assumptions from previous 
works [6—8]. These are: 

1. The fluid flow is laminar and fully developed. 

2. The thickness of the cylinder or the displacer is greater than the 
thermal penetration depth. 

3. Radiation heat transfer is neglected. 

4. Axial temperature gradients in the displacer, in the cylinder and 
in the gap are constant and equal to each other. This gradient is 
approximated by F = Th - TilLd- 

However, 

• Contrary to Baik [6] and Chang [7], the gas is not assumed 
incompressible and pressure gradient effect in the gap is taken 



























M.T. Mabrouk et al. / Energy 72 (2014) 135-144 


137 



into account but pressure variations in compression and 
expansion spaces are considered small compared to their mean 
value. This is generally the case for Stirling engines to avoid 
Hysteresis losses. 

• contrary to Kotsubo [8] the leakage flow rate is not considered 
known. 

• contrary to all of them we choose to use cylindrical coordinates 
instead of parallel plates geometrie. 

• contrary to all of them, we dont use a first order approximation 
to model fluctuations and this allows us to introduce our model 
into a realistic engine model. 


In a Stirling engine gas velocity and displacer velocity are fluc¬ 
tuating with zero mean value. x d (t) has also a zero mean value in 
our coordinate system (Fig. 2). However, gas pressure and different 
temperatures are composed of a time averaged and time- 
dependent parts. 

Respectively to previous assumptions, we have: 

• In the fluid gap 


p(z, t) = p 0 (z) + P(z,t) 

(2a) 

u(z, r, t) = u(z, r, t) 

(2b) 

T f (z,r, t) = T ref - Tz + T f (r, t) 

(2c) 

• In the displacer 

u d (t) = u d (t) 

(3a) 

x d (t)=x d (t) 

(3b) 

T d (Zd , r d , t ) = Tref - Tz d + T d (r d , t) 

= T ref -rz + rx d (t) + T d (r, t) 

(3c) 


• In the cylinder 

T c (z,r,t) =T ref - rz + f c (r,t) (4) 

where T re f is a reference temperature, which is simply selected as 
the cylinder temperature at a location with z = 0. 



Fig. 3. CPU-3 engine (0 type with Rhombic mechanism). 


2.3. Governing equations 


According to previews assumptions, momentum and energy 
equations for the flow in the gap are given in the literature 11 ]: 


dll 1 dp (d 2 u 1 du\ 

at pf 0 dz + Vf0 ^ar 2 + r dr J 


( 5 ) 


37> 37>o._ 

at az 


TfoPo dp 


“/ o 


(d 2 f f 1 dt f 
l ar 2 + T~dF 


( 6 ) 


P/oCp/o St 

In order to take into account heat conduction in the cylinder and 
in the displacer, we write governing equations in these parts: 


dT c 

~dt 

dTd 

at 


= “c o 


= “do 


a 2 t c 1 3T C \ 
dr 2 + r dr J 

'd 2 T d 1 dT d \ 
dr 2 + r dr J 


(7a) 

(7b) 


2.4. Boundary conditions 


We use non-slip boundary conditions for gas velocity: 


u(R d ,t ) = u d (t) 

(8a) 

u(R c , t) = 0 

(8b) 


According to assumption 3, the cylinder and the displacer can be 
considered as semi-infinite walls. Thus, temperatures at r = R d - S d 
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Table 1 

Geometric parameters of the GPU-3 engine. 


Parameter 

Value 

Clearance volumes 


Compression space 

28.68 cm 3 

Expansion space 

30.52 cm 3 

Swept volumes 


Compression space 

113.14 cm J 

Expansion space 

120.82 enr 

Heater 


Tube number 

40 

Tube inside diameter 

3.02 mm 

Tube length 

245.3 mm 

Void volume 

70.88 cm 2 

Cooler 


Tube number/cylinder 

312 

Tube inside diameter 

1.08 mm 

Length of the tube 

46.1 mm 

Void volume 

13.8 cm 3 

Regenerator 


Diameter 

22.6 mm 

Length 

22.6 mm 

Wire diameter 

40 pm 

Porosity 

0.697 

Unit number/cylinder 

8 

Thermal conductivity 

15 W/mK 

Void volume 

50.55 cm 3 

Displacer's length 

7.01 cm 

Cylinder's diameter 

6.96 cm 


and r = R c + S c are steady and not affected by the perturbation 
caused by the fluid flow. Temperatures at these radius are linear: 


T c (z, R c -f- <5 C , t) — T re f — rz 

(9a) 

Td( z d, R d - 5 d, t) = T ref - rz d 

(9b) 

Eqs. (9a) and (9b) allow us 
r = R d - S d and r = R c + S c : 

to write boundary conditions at 

Tc(R d + Sc, t) = 0 

(10a) 

Td(Rc - Sd, t) = 0 

(10b) 

Finally, we have a temperature and heat flux continuity condi¬ 
tions at r = R c and r = R t p 

T c (z,R c ,t) = T f (z,R c ,t) 

(Ha) 

T d (z,R d ,t) = T f (z,R d ,t) 

(lib) 

, 3 Uz,R d ) , 9T f (z,R d ) 

*d0 dr - *f0 dr 

(11c) 

dT f (z,Rc) . dT c (z,R c ) 
dr = Ac0 3r 

(lid) 


Table 2 

Operating conditions of the considered engine. 


Parameter Value 


Working gas 

Helium 

Frequency 

41.7 Hz 

Hot temperature 

977 I< 

Cold temperature 

288 1< 

Mean pressure 

4.13 M Pa 

Experimental power 

3.9 kW 


3. Analytic solutions for governing equations 


In order to solve governing equations in a general case we use 
Fourier series. In fact, all variables of the problem are periodic and 
can be expressed as Fourier series [12]. When taking complex no¬ 
tation, we can write: 

CO 

p\z, t) = p 0 {z)+ ^p n {z)e iut 

n= 1 

(12a) 

CO 

u{z,r,t) = ^u;(z,r)e iwt 

77=1 

(12b) 

CO 

u d = E U dn eiU 

77=1 

(12c) 

CO 

x d = Ev 1 "" 

77=1 

(12d) 

oo 

T} (z,r, t) = Tj 0 (z) + E T/h(r)e iwt 

n=l 

(12e) 

oo 

T*(z, r ,t) = Tco(z) +J2 r cn(r)e iMt 

n =1 

(12f) 

oo 

T d {z, r,t) = T d0 (z) + E Kn( r ) ei<ot 

(12g) 


n=t 


3.1. Velocity distribution of the gap fluid 

By reformulating Eq. (5) in complex notation and regrouping 
terms of the same order. We have for the nth order: 

d 2 K 1 du* n into ♦ _ 1 dp; 

Or 2 r dr vf 0 n p/o^fo dz 

Equation (13) is a modified Bessel equation. Then, it can be 
solved analytically. Its solution is [13]: 

1 On* 

= -s—=— ^ L + C„/ 0 ( e „r) + D R /o( en r) (14) 


a >0 - GHS SHI VI m 

me —► Qk *■ Q’ *“■ Qh me F 

_iJJ 



Fig. 4. Schematic model for five compartments of the Stirling engine and temperature 
distribution in each compartment as in Ref. [4]. 
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Table 3 

Set of differential and algebraic equations for the ideal Adiabatic Stirling engine model. 
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Type 

Equation 

Pressure 

dp = -yp(dV c /T ck + dV e /T he )/V c /T ck + y(V k /T k + V r /T r + V h /T„ + V e /T he ) 

Masses 

™k = pV k /RT k 
m r = pV r /RTr 
m/i = pV h /RT h 

m e = M - (m c + m k + m h + m e ) 

Temperatures 

T c = pVc/Rmc 

T e = pV e /Rm e 

Mass variations 

dm c = (pDV c + V c Dp/y)/RT ck 
dm k = m k dp/p 
dm r = m r dp/p 
dm h = m h dp/p 

Mass flow rates 

m ck = -dm c 
m kr = m ck - dm k 
m r /, = m kr - dm r 

m/ic = nirf, - dm e 

Conditional temperatures 

T ck = T c if m ck > 0 

T ck = T k if m ck < 0 

Ti,e = T h if m hs > 0 

The = T e if m /le < 0 

Energy 

dW = p(dV c + dV e ) 
dQk = V k dpc„/R - c p (T ck m ck - T kr m kr ) 
dQr = V r dpc t ,/R - c p (T kr m kr - T rll m rh ) 
dQh = V h dpc„/R - c p (T rh m rh - T he m he ) 


where /o et Jo are respectively modified Bessel functions of the first 


and second species and e n = Jino>/vf 0 . 

Constants C n and D„ are determined using non-slip boundary 
conditions: 



cranck angle <|> (°) 


Fig. 5. Chart flow of the simple analysis and introduction of the displacer gap losses 
model. 


Fig. 6. Position and velocity of the displacer. 
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to the displacer's velocity. Gas velocity fluctuation is the conse¬ 
quence of these two causes. The real velocity is derived from the 
complex one using the following relation: 

CO oo 

u(t) = ^K(u*)cos(n&jf) - ^3(u*)sin(ntot) (17) 

n=l n= 1 


3.2. Expression of leakage mass flow rate 


Leakage mass flow can be calculated by integrating the product 
of the density and the velocity over the area of the gap. In complex 
notation, it comes: 


Rc 

mj(t) =2np f0 J u(r,t)rdr (18) 

Rd 

We give an analytical expression of this complex mass flow rate: 


min = - 


8 z 
icj 


Rc 


2tcp/qD 

£n 


- R 2 d ) + 2 ^ p I° C l [R c/l (e n R c ) - R d f (e n R d )] 
-[RcJi(e n R c ) ~ R p h(e n R d )} 


(19) 


The real leakage mass flow rate is then given by the following 
expresion: 


rn/(t) = J 2 ^( lh \n)cos(ii 0 )t) - ^3(mjn)sin(nwt) (20) 

n=l n=l 


3.3. Temperature distribution 

By reformulating Eqs. (6), (7a) and (7b) in complex notations 
and regrouping terms of the same order n, we obtain: 


dT rn e 2 prT = zi ^Y n Pru n - e 2 n Pr^-p n 

r dr n f n n<o dz n n n Pf0 C pf0 Fn 


> 

dr 2 


(21a) 


a 2 t* 

u 1 cr 

dr 2 


1 37^ 
r 3r 


+ ~ Z" — t 'nR ra cT cn — 0 


(21c) 


where Pr is the gas Prandtl number, a c = ctfolu c o and a d = o/o/ado- 
Solutions of Eqs. (21a), (21b) and (21c) are: 


’9a, 


T fii = ~ 


PfoinroY 


*_ +i 

Pr - 1 noj 


R dn — Rnh 0 ^nx/PtCdtj + LnJo (rn'/Rta d r S J 
T cn = G n Io(e n \/PE<7 C f) + HnJo f e ny/P ra c r 


1 Pr - 1 nw 

■r) 

TfofioPn 


PfoCpfo 



(22a) 

■) 

(22b) 

■) 

(22c) 


where E n , F„, G n , H„, K n , and L n are integration constants which can 
be calculated using boundary conditions given by Equations (9a) 
and (9b)—(lid). Writing these equations using Fourier series we 
have a system of equations for each order n: 


T dn( R d ~ 5 d) = 0 

(23a) 

dT dY R d) *f0 dT fn( R d ) 

0r A d0 dr 

(23b) 

^dn(Rd) + Ix dn = Rfn( R d) 

(23c) 

dT* cn (Rc) */0 3 T fr( R c) 

dr A c0 dr 

(23d) 

= Tfn(Rc) 

(23e) 

Ren (Re + <5 C ) = 0 

(23f) 


This equations system is linear and can be solved by hand, but 
solutions obtained are lengthy. Thus, we choose to solve it 
numerically. 


a 2 r 


dn 


i ar; 


dr 2 


H— ■ 
r 


dn 

dr 


- epr'dTdn = 0 


(21b) 3.4. Calculating displacer gap losses 



At an arbitrary axial position z, the instantaneous local enthalpy 
flow rate from the high to the low temperature side in the gap and 
the displacer are given by Ref. [8]: 

• In the gap 

R c 

H g (t) = 2t i J ( p f0 C pf0 T f u + (l - T f0 f3 0 'jpujrdr (24) 

R„ 

• In the displacer 

Rd 

H d (t) = 2n J (p d0 C pd0 T d u d ydr (25) 

Rd-dd 


Fig. 7. Pressure drops in the heater, the regenerator and the couler and the total Notice that real gas effects are taken into account by the term 

pressure drop p c - p e . (1 - Tfo(3o)pu which is null for ideal gas. 
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Fig. 8. Gas velocity distribution at different crank angles. 


W) = ^ j H g (t)dt+ 2 n j H d (t)dt (26) 

o o 


4. Case study 

In order to calculate leakage losses in a real Stirling engine, the 
GPU-3 engine [14] is chosen as a case study. This /J type Stirling 
engine was developed by "General Motors" in 1965. It is shown in 
Fig. 3 and its specifications and operating conditions are gathered 
in Tables 1 and 2. Prior to evaluate its displacer gap losses, we have 
to develop a thermal model for the engine. 

4.3. Thermal model of the engine 



Fig. 9. Leakage mass flow through the displacer's gap. 


Displacer gap losses are the sum of "enthalpy pumping" through 
the displacer's clearance and "shuttle heat transfer". "The enthalpy 
pumping" is defined as the cycle averaged net enthalpy flow rate 
from the high to the low temperature side in the fluid gap. The 
"shuttle heat transfer" is defined as the cycle averaged net enthalpy 
flow rate in the displacer. Then, displacer gap losses are given by 
Refs. [6-8]: 


The simple analysis is used in order to model the engine. This 
analysis is a decoupled losses approach based on the ideal adiabatic 
model proposed by Urieli and Brechowitz [4[. Stirling engine is 
divided, as illustrated in Fig. 4, on five compartments: an expansion 
space, a heater, a regenerator, a cooler and a compression space. 
There are four interfaces between them through which the working 
gas passes from a compartment to another. Governing equations of 
the adiabatic model are obtained by applying energy and mass 
balances and the equation of state on each compartment. The set of 
equations of the ideal adiabatic model are gathered in Table 3. It 
contains non-linear differential equations to be integrated numer¬ 
ically. We use a classical fourth order Runge-Kutta method for the 


3.56 
3.54 
3.52 
"g" 3.5 


“ 3.48 

t5 

ro 

CC 3.46 
3.44 
3.42 
3.4 



o Displacer ((f>=60°) 
o Gap (<J>=60°) 
o Cylinder (<|)=60 o ) 

□ Displacer ((f)=120°) 

□ Gap (<f>=120°) 

□ Cylinder (t|>=120°) 

* Displacer ( 4>=1800°) 

* Gap (<f>=180°) 

* Cylinder (<(>=180°) 

* Displacer (((>=240°) 

* Gap ((j)=240 o ) 

* Cylinder ((|>=240 o ) 

* Displacer (<t>=300°) 

* Gap (<t>=300°) 

* Cylinder ((j>=300°) 


500 550 600 650 700 750 800 

Temperature (K) 



450 500 550 600 650 700 750 800 

Temperature (K) 


O Displacer (<J>=60°) 
O Gap (<|>=60 o ) 

O Cylinder (t))=60 o ) 

□ Displacer (<(>=120°) 

□ Gap (<t>=120°) 

a Cylinder ((j)=120°) 

* Displacer (<|>=180°) 

* Gap (<t>=180°) 

* Cylinder ((j)=180°) 

* Displacer (t()=240 o ) 

* Gap (<j)=240°) 

* Cylinder (t))=240 o ) 

* Displacer (<|)=300 o ) 

* Gap (<|>=300 o ) 

| » Cylinder ((|)=300 o ) 


500 550 600 650 700 750 800 

Temperature (K) 



Fig. 10. Radial distribution of gas temperature in the gap at z = 0 and at different crank Fig. 11. Radial distribution of the temperature in the displacer, in the gap and in the 
angles. cylinder for z = 0 and at different crank angles and for two clearance thicknesses. 
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integration starting from an initial condition defined using the 
Schmidt analysis [2], 

In the simple analysis, heat transfer, pressure drop and effec¬ 
tiveness of the regenerator are taken into account following the 
scheme presented in Fig. 5. 

Input data for the simple analysis are geometric parameters and 
operating conditions. Instantaneous positions of the piston and the 
displacer are deduced from a cinematic study of the engine’s drive 
[4], Instantaneous position and velocity of the displacer are given 

in Fig. 6. 

The pressure difference between the compression and expan¬ 
sion spaces is the sum of pressure drops in the heater, the regen¬ 
erator and the cooler with respect to the direction of the flow Eq. 

(27): 

Pc(t) ~ Pe(t) = J2 Sign(U;)Apj (27) 

ie[h,r,k] 


where i = h for the heater, i = k for the cooler and i = r for the 
regenerator. u, is the gas velocity and Ap,- is the pressure drop 
calculated by: 


A Pi = 


fiLiPiuj 

2D, 


(28) 


fi is the friction factor, D ,- is the hydraulic diameter, L, is the 
length and p,- is the gas density. 

For the hot and the cold heat exchangers, the friction factor is 
equal to 64/Re in laminar regions and is calculated using a corre¬ 
lation given by Petukhov [15] in turbulent regions. For the regen¬ 
erator, the friction factor is calculated using a correlation developed 
by Tanaka [16]: 


• for the heater and the cooler: 


64/Re if Re,- <2300 

(29) 

(0.79 In (Re,) - 1.64)'~ 2 if Re,- > 2300 

(30) 


• for the regenerator: 


fr 


175 

RSr, max 


+ 1.6 


(32) 


These pressure drops and total pressure drop (p c (f)-Pe(f)) are 
given in Fig. 7 where numerical data are fitted using third order 



0 50 100 150 200 250 300 350 

Cranck angle 4> (°) 



Fig. 13. Enthalpy flow rate in the displacer at z = 0 and its mean value (shuttle heat 
transfer). 


Fourier series. Please notice here that the subscript c stands for 
"compression space". 


4.2. Velocity and temperature distributions and displacer gap losses 

We adopt arbitraly z = 0 and we define the pressure's gradient at 
this location by a spatially averaged value: 3p(t)/9z = (p c (t) - Pe(t))/ 
La. For the instantaneous pressure we define it as the arithmetic 
average between pressures in the compression and expansion 
spaces p(t) = (p c (t) + p e (t))/2. This is possible without a great error 
because pressure's spacial variation in the gap is very small 
compared to its temporal variation. 

The instantaneous velocity distribution in the gap between the 
displacer and the cylinder at this location is then calculated using 
Eq. (17). In Fig. 8, we show this distribution at different crank angles 
for a clearance thickness of 50 microns. Notice that it respects 
boundary conditions. Leakage mass flow rate is then calculated 
using Equation (20) and shown in Fig. 9 for two clearance thick¬ 
nesses (50 and 80 microns). 

After solving the system of equations constituted by 
Eqs.(23a)— (23f) for each order, we can give the temperature dis¬ 
tributions in the displacer, in the fluid gap and in the cylinder. Gas 
temperature distribution in the gap is given in Fig. 10 for a clearance 
thickness of 80 microns and at different crank angles. In Fig. 11 we 
can find temperature distributions in the three regions (cylinder, 



Fig. 12. Enthalpy flow rate through the gap: Enthalpy without real gas effect, real gas Fig. 14. Effect of the clearance thickness on the shuttle heat transfer and the enthalpy 
effect and total enthalpy flow rate. pumping through the displacer's gap. 
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Fig. 15. Variation of the enthalpy pumping, the shuttle heat transfer and the total 
displacer gap losses as a function of the frequency. 

displacer and fluid gap) at different crank angles and for two 
different clearance thicknesses (50 microns and 80 microns). One 
can notice in this figure that temperature continuity is respected. 

Enthalpy flows in the gap and in the displacer at z = 0 are 
calculated for a clearance thickness of 50 microns using Eqs. (24) 
and (25). Fig. 12 shows the enthalpy flow through the gap 
without real gas effect (the term p 0 C p oTju only), the effect of real gas 
(the term (1 - Tjo/3o)pu ) and their sum which represents total 
enthalpy pumping through the gap. We note that for this particular 
case cycle averaged real gas effect represents around 17% of total 
cycle averaged enthalpy flow pumped through the clearance 
(31.1 W/180.1 W). 

Fig. 13 shows the enthalpy flow in the displacer at the same fixed 
axial coordinate z = 0 and the cycle averaged value of this enthalpy 
flow is equal to 223 W which give us total displacer gap losses of 
404.5 W. 

4.3. Effect of the clearance thickness 

Enthalpy pumping and shuttle heat transfer was investigated by 
varying the clearance thickness and results are given in Fig. 14. 
Shuttle heat transfer from our model was compared to analytical 
solutions given by Baik [6] and Chang [7], We notice that for small 
clearance thicknesses where the effect of the enthalpy pumping is 
small, results from Baik's formula is very close to ours. 

However, his results become different from ours when the 
clearance thickness increases. Our results show that shuttle heat 
transfer decreases more rapidly than predicted by both Baik and 
Chang and may even be negative. This is due to the effect of the 
leakage gas flow in the gap wich was not taken into account by the 
former authors. In fact, the displacer exchanges heat with the fluid 
and returns back a part of the energy pumped through the gap. This 
makes sense since we notice that gas velocity is always opposite to 
displacer velocity (Fig. 8) and this is due to the effect of the pressure 
gradient in the gap that we take into account. So, when clearance 
thickness rises, enthalpy pumping rises and shuttle heat transfer 
falls. The consequence of this opposite behavior is that total dis¬ 
placer gap losses present a minimum situated in the studied case 
between 50 and 70 microns. 

4.4. Effect of the engine's speed 

It's known that the pressure drop in Stirling engines is more 
important at high speeds. So, one can expect that engine's speed 
affects the enthalpy pumping. This fact is confirmed by Fig. 15 
where we draw enthalpy pumping as a function of the frequency. 


Fig. 15 shows that the enthalpy pumping rises as the frequency 
rises. 

For the shuttle heat transfer, we notice that its behavior is 
notably different from Baik’s and Chang's solutions which predicts 
that shuttle heat transfer increases slightly up to 15 Hz and is 
constant for higher frequencies. Our results (Fig. 15) show that 
shuttle heat transfer decreases when frequency increases. This is 
due once again to the effect of the leakage flow in the gap which 
become more important at high speeds. 

Consequently, total displacer gap losses present a minimum 
situated in the studied case between 25 and 35 Hz. 

5. Conclusion 

This study is devoted to "Displacer gap losses" in (j and y Stirling 
engines. "Displacer gap losses" are the sum of "enthalpy pumping" 
and the "shuttle heat transfer". Enthalpy pumping is due to the 
leakage flow between compression and expansion spaces. The 
shuttle heat transfer is due to the displacer motion between hot 
and cold regions. 

In this work, we develop, in the first part a model for the pul¬ 
sating fluid flow through the displacer clearance where we take 
into account the pressure difference between compression and 
expansion spaces. Gas velocity and temperature distribution in the 
gap are derived from this analytical model. Conduction heat 
transfer equations in the cylinder and in the displacer are solved 
respect to heat and temperature continuity boundary conditions 
with the gas flowing in the gap. 

In the second part, we introduce our model in a global Stirling 
engine model. GPU-3 is chosen as a study case. In respect of the 
crank angle, we give gas velocity and temperature distributions in 
the displacer, in the gap and in the cylinder. Then, we investigate 
the behavior of the gap losses by varying the clearance thickness 
and engine's frequency. 

For clearance thickness, our results demonstrate that displacer 
gap losses show a minimum value. For our study case, this mini¬ 
mum is located in the range between 50 and 70 microns. This is a 
new result compared to former studies. 

For frequency, our results demonstrate that displacer gap losses 
show a minimum value. For our study case, this minimum is located 
in the range between 25 and 35 Hz. This is also a new result 
compared to existing studies. 
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